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Study of the hadronic matrix elements can provide not only tests of the QCD sector of the Stan- 
dard Model (in comparing with existing experiments) but also reliable low-energy hadronic quan- 
tities applicable to a wide range of beyond-the-Standard Model scenarios where experiments or 
theoretical calculations are limited or difficult. On the QCD side, progress has been made in the 
notoriously difficult problem of addressing gluonic structure inside the nucleon, reaching higher- 
Qh Q 2 region of the form factors, and providing a complete picture of the proton spin. However, 

even further study and improvement of systematic uncertainties are needed. There are also pro- 
posed calculations of higher-order operators in the neutron electric dipole moment Lagrangian, 
which would be useful when combined with effective theory to probe BSM. Lattice isovector 
tensor and scalar charges can be combined with upcoming neutron beta-decay measurements of 
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the Fierz interference term and neutrino asymmetry parameter to probe new interactions in the ef- 
fective theory, revealing the scale of potential new TeV particles. Finally, I revisit the systematic 



uncertainties in recent calculations of gA and review prospects for future calculations. 



The XXX International Symposium on Lattice Field Theory 
June 24-June 29, 2012 
Cairns, Australia 



* Speaker. 
1"NT@UW-12-19 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http : / /pos . sissa . it / 



Lattice Hadron Structure: Applications within and beyond QCD 



Huey-Wen Lin 



1. Introduction 

Calculation of the low-energy hadronic matrix elements is important not only to make Standard 
Model predictions for comparison with experimental values but also to provide the valuable inputs 
needed to assist in the searches for new physics. Flavor physics from lattice QCD has been playing 
an important role for the Standard Model inputs in heavy meson sector; B- and D-meson decay 
constants, B-mesons mixing parameters and more. For lighter mesons, such as pion and kaon 
systems, lattice QCD can determine many related quantities pretty reliably. However, when it 
comes to the nucleon, difficulties arise; this is frustrating, since unlike the mesons, which exist 
only ephemerally, nucleons are one of the building-block particles that compose us, and most 
measurements to understand Standard Model and probe new physics will involve interactions with 
nucleons. Given the difficulties and importance of nucleon quantities, in this review, I will focus 
on the nucleon, since there are still many challenges that need to be resolved and exciting new 
developments. 

Nucleon matrix-element calculations on the lattice are harder in general because of the fol- 
lowing problems. Firstly, we face a signal-to-noise problem. In Euclidean space the nucleon 
signal-to-noise ratio scales like exp(— M# + 3m K /2) with Euclidean time; we need to extract the 
ground-state signal before noise becomes overwhelming. The first excited-state 1440) is nearby; 
therefore, we have to be careful to distinguish the ground and excited signals to avoid contamina- 
tion and to wait for larger time where the excited-state signal is relatively diminished, if necessary. 
However, this is at odds with the first problem; as we wait for larger time to avoid excited states, 
the signal-to-noise ratio decreases. Thus, we need high statistics and longer trajectories to have 
reliable numbers for experiments. Second, chiral perturbation theory (ChPT) is more difficult, be- 
cause the A resonance is nearby. As a result, there are multiple expansions for ChPT, and some of 
them may have poor convergence, bringing the validity of chiral forms into question at the heaviest 
relevant pion mass. Thankfully, this issue may soon become moot in the physical pion-mass era. 
Also, larger volumes are required to avoid systematic uncertainties due to the finite volume. Low- 
statistics and heavy-pion studies of finite-volume effects cannot be directly applied to the case of 
light pion mass. And further, gauge ensembles are not always generated with nucleons in mind, 
resulting in suboptimally small volumes and low numbers of configurations. 

Given the limited space, I will concentrate on reviewing a selection (colored by personal in- 
clination) of exciting topics and progress since the last lattice conference plenary talk on the same 
topic. I start with hadron structure, focusing on new calculations of the gluonic momentum fraction 
and calculations of how pion structure is modified in a pion medium. Then, I will give two examples 
of calculations with applications to searches for physics beyond the Standard Model (BSM): the 
neutron electric dipole moment and new tensor and scalar interactions in neutron /3 decay. Lastly, I 
will review recent calculations of the nucleon axial charge and carefully consider its systematics. I 
will present more details on excited-state contamination, and I conclude that a better evaluation of 
the finite-volume effects is urgently needed to understand the current discrepancy between lattice 
calculations and the experimental value. 

2. Probing Inside Hadron 

The picture inside the hadron is usually described by studying longitudinal quark distributions 
in momentum space (structure functions, such as quark and gluon momentum fractions (x n ) q , g 
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to the n moment), or studying transverse quark distribution in coordinate space (through form 
factors). However, the ambitious generalized parton distributions (GPDs) give a fully-correlated 
quark distribution in both coordinate and momentum space. Given the limited space, only some 
selected interesting updates will be presented here. Please also refer to the previous lattice hadron 
structure reviews [1, 2, 3, 4, 5] for omitted definitions. 

2.1 Structure functions 

One of the fundamental questions in QCD is by what proportion the gluons and quarks make 
up the hadron and give their masses. The mass of pion and nucleon differ by an enormous amount 
despite there being only an increase of one quark in the hadron; the gluons must play a significant 
role. One way to directly study the question is to look at the gluon momentum fraction (x) g (and 
its higher moments). Due to rotational symmetry breaking on the lattice, we can use either fff = 
tr(<f x St) or G\ = § tr(<? 2 + S 2 ) for (x) g However, gluonic structure has been notoriously difficult 
to calculate with reasonable signals in lattice QCD, even for just the first moment. Despite these 
difficulties, the gluonic structure can be studied in the quenched approximation, and these new 
updates provide approaches and successful demonstrations that give some hope that the problem 
can be addressed. Both #QCD [6] and QCDSF [7] have made breakthroughs with updated studies 
of gluonic moments in the nucleon using Wilson gauge actions (j3 = 6.0, aT x fa 2 GeV) with Wilson 
(m K £ {480,650}) and nonperturbative 0(a) -improved clover {m n G {480,650}) fermion actions, 
respectively. 

#QCD [6] uses a massless overlap Dirac operator for the gauge-field tensor as ti s (a^ v D 0V ) oc 
a 2 fuv to compose the operator G s = F^F^. The calculation is done using 500 gauge configu- 
rations, and two Z4 noise sources to stochastically estimate D ov for all space points at a given time. 
The same technique can be applied to look for general gluonic matrix-element contributions. 

QCDSF uses the Feynman-Hellmann theorem by adding 



into the total QCD action. Then they calculate the nucleon mass (using standard two-point cor- 
relators, as in the nucleon sigma term) as a function of A, which couples to the relevant gluonic 
operators; the derivative of the mass with respect to A contains the wanted gluonic matrix element: 

M - -2 J_ dM "W 
Wg ~ 3M N dX A=0 " 

Both works use similar statistics, 0(500) configurations, and result in similar signal-to-noise 
ratio of the final result when extrapolated to the physical pion mass, about 15%. The QCDSF 
approach is interesting and computational cheap, but requires new configurations for each operator. 
One may need reliable reweighting techniques when moving to dynamical lattices. The #QCD 
approach is more robust; the same improved field operator can be used for any gluonic operator. 
However, it requires all-to-all overlap Dirac operators at any given time, so it is relatively expensive. 
However, one may be able to reuse D ov for quark disconnected contributions and general gluonic 
operators. 

Another interesting development is the study of how the hadron structure changes with in the 
presence of a medium. One of the most exciting topics in nuclear physics is the study of how 
hadronic properties are modified by the nuclear environment and how such modifications affect 
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Figure 1: (left) The pion quark momentum fraction in n + , expressed as a ratio with its value in the vacuum, 
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I (x) 71 ' as a function of pion density, (right) Total, gluon and connected up/down quark angular mo- 
mentum to proton spin from jQCD collaboration. Disconnected up/down and strange quark contributions 
are shown in the bottom panel. 



the properties of nuclei. In 1983, the European Muon Collaboration (EMC) reported surprising 
deep inelastic scattering (DIS) measurements of the structure function F2 as a function of Bjorken 
x; the measurement changed significantly between heavy nuclei and deuterium. Since then, many 
experiments have refined this measurement and looked at a variety of other isotopes and light 
nuclei. Many models were devised to try to explain the effect, and although many suggested 
factors appear to play some role, there is no universal understanding of the EMC effect yet. Such 
understanding would also be important to provide a better picture of how in-medium properties 
change to provide a baseline to constrain BSM physics, for example, in the NuTeV anomaly weak 
mixing angle experiment. This measurement lies 3 sigma away from the SM expectation, but the 
discrepancy may be explainable within QCD by environmental nucleon-effect corrections. 

Unfortunately, multibaryon systems are complicated to calculate in LQCD; the noise increases, 
similarly to the sign problem. Therefore, we start with pion structure in a pion (7T + ) medium. 
Multipion systems suffer from significant thermal effects due to the finite time extent of the box, 
due to pions going "around the world", and the problem worsens for light ensembles. To extract 
the energy of the multipion state, we need to not only consider the contribution of the n-% system, 
but also all combinations of lower numbers of pions that can travel forward and backward ("around 
the world") through the periodic time extent of the lattice volume. 

This first lattice-QCD attempt to measure EMC effects for the pion momentum fraction in a 
pion medium uses pion masses ranging 290-490 MeV at 2 lattice spacings [8]. On each ensemble, 
3-5 source-sink separations are used to take the thermal-state degrees of freedom into account. The 
updated preliminary results for the ratio of the quark momentum fraction of the pion in medium 
versus the vacuum (x) n ' N / (x) n ^ are shown on the left-hand side of Fig. 1 as function of pion 
density. The gray band is the preliminary extrapolated physical limit, indicating strong medium 
corrections to the first moment of the pion quark momentum fraction. 

2.2 Form Factors 

The low-<2 2 region of hadron form factors is a main focus of lattice calculations. We should 
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devote more effort and resources to improving the lattice-QCD calculation of the charge radii; 
currently, the majority of calculations of the charge radii are about a factor of two smaller than the 
experimental value. Given the current tension and on-going theoretical and experimental efforts 
devoted to this quantity, we should seek better constraints on the radii. 

Another interesting quantity is the induced pseudoscalar coupling constant g* p , which is impor- 
tant for muon-related physics; many of these experiments are designed to search for new physics. 
g* P comes out as side product of the nucleon axial-current form factor, but in experiment, the in- 
duced form factor is defined at transfer momentum 0.88m 2 . Currently only RBC/UKQCD are 
reporting this number, and the constraints are poor at the moment. (See Fig. 5 in Ref. [9] for a 
summary plot.) Recent experimental results from MuCap [10] have updated the g* P precision to 
less than 1%, 8.06(55). 

Studying the high-<2 2 region helps us investigate the distribution of charge density within the 
hadron. Even though at sufficiently large momenta, perturbative QCD should be able to describe 
the physics well, recent pion-to-two-photon data fail to converge to the asymptotic limit even at 
40 GeV 2 . An independent reliable nonperturbative QCD method would help. For example, work 
in Ref. [11] extends the available Q 2 region by keeping nucleon operators that overlap both with 
zero-momentum and highly boosted ground-state nucleons. This helps out with the signal-to-noise 
at a given lattice spacing. Taking the large-<2 2 form-factor data from Ref. [11], we can map out 
the transverse charge-density distribution as a function of the impact distance b in a polarized 
nucleon [12, 13, 14]: 

Pr(b) = [^MbQ)Fi(Q 2 ) + sm^)^^MbQ)F 2 (Q 2 ), (2.1) 

where 7o,i are Bessel functions. We can perform this integral numerically, using the lattice F\^{Q 2 ) 
obtained by extrapolating the fit form to the physical pion mass. Figure 2 shows the results for the 
proton and neutron in the two-dimensional impact plane using lattice inputs. There are positive 
and negative charges surrounding the center, which in the neutron sum to zero. Ref. [11] also show 
that using Q 2 < 3 GeV 2 will lead to significant deviations from the true short-distance behavior; 
compare with, for example, 6 GeV 2 . Note that the asymmetry in the distribution for a polarized 
nucleon is due to the relativistic effect of boosting the magnetic moment of the baryon. This induces 
an electric dipole moment that shifts the charge distribution. 

2.3 Proton Spin 

One of the fundamental questions in QCD is how the proton's spin 1/2 is distributed among 
its constituents. The most naive intuition is that the three quarks carry the spin. Early experiments 
found the quark contribution makes up less than half of the total spin; this contradiction was dubbed 
the spin crisis, and many experiments and models attempted to address it. Now we better under- 
stand that significant contributions come from both intrinsic spin and orbital angular momentum. 
Ji [15] gave a set of covariant operators for these that we can calculate on the lattice. QCDSF, 
LHPC and ETMC have been reported quark connected pieces for spin and orbital angular momen- 
tum for many years now (see reviews in Refs. [1, 2, 3, 4, 5, 16]). These results are renormalized 
in the MS scheme at a scale of 2 GeV and shown on the right-hand side of Fig. 1 . The extracted 
values are the stars on the left. As it turns out, the quarks contribute less than 50% of the total 
nucleon spin; the total quark orbital angular momentum contribution is consistent with zero, so we 
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Figure 2: The transverse charge densities in a polarized proton (left) and neutron (right) with impact pa- 
rameter b x y in fm. The color black indicates near-zero values, and purple, orange and white are increasingly 
positive (figures now in Exploring the Heart of Matter in the National Academies Press). 

must conclude that a majority of the proton spin comes from gluons, which is quite a surprise. The 
goal is still moving to lighter pion masses, addressing disconnected diagrams and evaluating gluon 
contributions that currently must be estimated from sum rules. 

%QCD report the first lattice-QCD calculation including a complete set of the disconnected 
contributions and gluonic operators (see Ref. [6] and private communication with Y. Yang). On 
the right-hand side of Figure 1, we see their result as a function of pion mass squared. The proton 
spin breakdown is: quark spin contributed 50(2)%, quark orbital angular momentum contributed 
25(12)% (mostly from the disconnected diagram, since they also see total zero quark contribution 
from the connected quark orbital angular momentum), and the total gluon spin contribution is 
25(8)%. However, the current calculation is still quenched Wilson action, and an update with 
overlap on DWF sea is in progress. It would be also interesting to see corroborating work from 
other collaborations using complete sets of the diagrams. 

3. Applications beyond QCD 

There are many opportunities to probe physics beyond the Standard Model (BSM) with lattice- 
QCD calculations of hadronic matrix elements. For example, we can calculate the leading hadronic 
vacuum contribution to the muon anomalous magnetic moment or evaluate processes involved in 
light-by-light contributions to it (by direct calculation, four-point correlators or meson-to-two pho- 
ton form factors). The current discrepancy between the world-average theoretical and experimental 
values is around 3^1 a; future improvement to the experimental data (reducing the error by a factor 
of 4 by 2015) will require further understanding and error-reduction on the nonperturbative QCD 
theoretical end. (See T. Blum in these proceedings for a detailed review.) The nucleon scalar den- 
sity, including its difficult disconnected piece, is an important low-energy constant needed as input 
for many BSM models of dark matter, such as supersymmetry (SuSy); R. Young in these proceed- 
ings will also cover it. In this section, I will consider examples such as the neutron electric dipole 



6 



Lattice Hadron Structure: Applications within and beyond QCD 



Huey-Wen Lin 



moment (nEDM) and non-V— A interactions in neutron beta decay. Here I will briefly mention 
the progress made in lattice-QCD nEDM calculations and lattice-QCD calculations relevant to the 
detection of novel interactions in neutron beta decay. 

3.1 Neutron Electric Dipole Moment (nEDM) 

The neutron electric dipole moment is a measure of the distribution of positive and nega- 
tive charge inside the neutron. To generate a finite nEDM, one needs processes that violate CP- 
symmetry, for example by adding a CP-odd 0-term to the Lagrangian. The SM value of this 
quantity is very small, 10~ 30 e • cm. Although experiments do not have the necessary precision to 
measure the SM value; many BSM models which predict values that are higher than the experimen- 
tal upper bounds have been ruled out. This includes some parts of the parameter space for certain 
SuSy models. The most common type of lattice-QCD calculation mainly focuses on the leading 6 
term, the CP-odd contribution; a current combined analysis gives 0(30%) in statistical error alone 
(see Fig. 1 in Ref. [16] for a summary plot). The precision of the lattice-QCD calculation needs to 
be greatly improved; all-mode averaging (AMA) has been proposed to highly improve the current 
statistics even at near-physical pion mass (see Ref. [17] and E. Shintani in these proceedings). It 
may also be possible to combine the nucleon matrix elements with higher-order operators, such as 
quark electric dipole moment (qEDM) and chromoelectric dipole moment (CEDM) operators, (see 
T. Bhattacharya, these proceedings) using effective field theory and proposed LQCD calculations. 

3.2 Non-V— A Interactions in Neutron Beta Decay 

The measurement of non-Standard Model contributions to precision neutron (nuclear) beta- 
decay measurements would give hints of potential BSM particles at the TeV scale; if new particles 
exist, their fundamental high-scale interactions could appear at low energy in the neutron beta- 
decay Hamiltonian as new terms: 

H en = G F ^ ept A x/n k + I^ SM 0! ept x0r rk ^) , (3-D 

where Gp is the Fermi constant, Jy-A is the left-handed current of the indicated particle, and the 
sum includes operators with novel chiral structure. The new operators will enter with the coeffi- 
cients £ that are related to the TeV scale of the particles, similar to how dimensionful Fermi constant 
gives hints to the masses of the W and Z bosons of the electroweak theory before they were discov- 
ered. In the context of the nucleon, the new operators of this Hamiltonian will introduce low-energy 
coupling constants: gj = (n\uGn V d\p), gs = (n\ud\p), and lattice QCD is a perfect theoretical tool 
to determine these constants precisely. The search for BSM physics proceeds experimentally by 
measuring the Fierz interference term and neutrino asymmetry parameter of neutron differential 
decay rate. Upcoming experiments at LANL (UCNB and UCNb) and ORNL (Nab) plan to mea- 
sure Obsm to the 10~ 3 level by 2013 (and maybe after a few years improve to 10~ 4 ). This only 
requires g$j to be known to 10-20% (after summing all the systematic uncertainties) to be useful 
as a BSM probe. For more details about experimental and theoretical work on this subject, refer to 
Ref. [9] and R. Gupta in these proceedings. 

Ref. [16] contains a summary of previous knowledge of these charges, experiments, models 
and previous lattice-QCD inputs, as applicable. Here we will concentrate on the updates available 
this year. PNDME (by Gupta) and LHPC (by J. Green) reported updated calculations of gs t r after 
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Figure 3: Global analysis of all Nf = 2+ 1(+1) lattice calculations of gj (upper-left) and g$ (upper- 
right) with m K L > 4 and m K T > 7 cut to avoid systematics due to small spatial or temporal extent. The 
leftmost points are the extrapolated values at the physical pion mass. The two bands show extrapolations 
with different upper pion-mass cuts: m\ < 0.4 and m| < 0.2. The m % L < 4 data points are marked faded 
within each calculation; the lattice-spacings for each point are denoted by a solid line for a < 0.06 fm, 
dashed: 0.06 < a < 0.09 fm, dot-dashed: 0.09 < a < 0.12 fm and dotted: a > 0.12 fm. (Right) e s -e T 
allowed parameter region using different experimental and theoretical inputs. The green and purple dashed 
lines are the constraints from first LHC run bound and near-term expectations with running to 2 GeV to 
compare with low-energy experiment. 



Ref. [9]. We use the chiral formulation given in Refs. [18] and [19] for tensor and scalar charges, 
respectively, to extrapolate to the physical pion mass (see the left-hand side of Fig. 3). We see that 
the PNDME points greatly constrain the uncertainty due to chiral extrapolation in both cases and 
obtain gy QCD = 0.978(35) and ^ QCD = 0.796(79) (statistical errors for now). The LQCD values 
are better determined than other theoretical estimations from various models, which give rather 
loose bounds on these quantities; for example, g$ is estimated to be between 0.25 and 1. 

These low-energy matrix elements, the tensor and scalar charges, can be combined with ex- 
perimental data to determine the allowed region of parameter space for scalar and tensor BSM 
couplings (denoted e). Using the g$j from the model estimations and combining with the existing 
nuclear experimental data 1 , we get the constraints shown as the outermost band on the left-hand 
side of Fig. 3. Combining the expected experimental data in 2013 and existing measurements, and 
again, using the model inputs of g$j, we see the uncertainties in £^7- are significantly improved. 
Finally, using our present lattice-QCD values of the scalar and tensor charges, combined with the 
expected 2013 precision of experimental bounds on deviation of these neutron-decay parameters 
from their SM values, we found the constraints on Ssj are further improved, shown as the inner- 
most region. These upper bounds on the effective couplings Esj will correspond to lower bounds 
for the scales Asj at 5.6 and 10 TeV, respectively, for new physics in these channels. 

'For example, nuclear beta decay + — > + transitions and other processes, such as j5 asymmetry in Gamow-Teller 
60 Co, longitudinal polarization ratio between Fermi and Gamow-Teller transitions in 114 In, positron polarization in 
polarized 107 In and beta-neutrino correlation parameters in nuclear transitions. 
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Figure 4: (Left) Collected experimental values used in PDG 2012 average (the band) and the latest UCNA 
(2012) measurements on gA', there has been a slow increase in gA/gv m the past 15 years. The lower panel 
shows the physical-pion extrapolated gA numbers collected from dynamical 2+1-flavor and 2-fiavor lattice 
calculations using (3(a)-improved fermions [20, 21, 22, 23, 24, 25, 26, 27, 28]. Note the change in scale 
and that most of the errorbars here are statistical only; a few calculations quoted systematics, which are 
added to statistical ones as outer errorbar bands, marked with dashed lines, (right) Published and ongoing 
gA calculations using 0(a)-improved dynamical fermions, plotted as a function of nv\. The lattice spacings 
are grouped as described in the captions of Fig. 3. LHPC clover m K T < 7 data points are excluded here. 



How do the constraints from high-energy experiments compare? As demonstrated in Ref. [9], 
neither CDF nor DO is sufficient to provide useful constraints. We can estimate the £5 7; constraints 
from LHC current bounds and near-term expectations through effective Lagrangian 

% = -^Vudiud^ePLVe) - ^jVudiuO^d^ea^Pjye), (3.2) 
A s A T 

where T)jj = ±1. The high-energy bound are scaled down to 2 GeV to compare with low-energy 
prediction; more details can be found in Ref. [9]. By looking at events with high transverse mass 
from CMS and ATLAS in the ev + X channel and comparing with the SM W background, we 
estimated 90%-C.L. constraints on 7]s.t based on the first LHC run, yfo = 1 TeV L = 10 ftT 1 (the 
green line) and for a near-future run yfs = 8 TeV L = 25 fb _1 (the purple dashed line) on the 
right-hand side of Fig. 3. 

Furthermore, for detection experiments that have not yet observed any events such as n-h os- 
cillations or proton decay, lattice-QCD calculations can provide low-energy constants to constrain 
the experimental search ranges. The potential to search for new physics using these precision nu- 
cleon matrix elements during the LHC era and in advance of Project-X experiments at FNAL make 
lattice-QCD calculation of nucleon structure particularly timely and important. 

4. Nucleon Axial Charge gA 

The nucleon axial charge gA is a very important fundamental measure of nucleon structure, 
with applications such as dictating the rate of proton-proton fusion, which is the first step in the 
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thermonuclear reaction chains that power low-mass hydrogen-burning stars like the Sun, being an 
important input to measurements of BSM physics, affecting the extraction of V u d and neutrinoless 
double-beta decay. Presently, gA is best known from the experimental measurement of neutron 
beta decay using polarized ultracold neutrons by UCNA collaboration, which dominates the PDG 
average with uncertainty at the 0.2% level [29]. (See the upper panel of the left-hand side of Fig- 
ure 4 for the collected experimental gA measurements used in PDG2012 and the recently updated 
UCNA number [30].) This uses some assumptions about new few-TeV particles that can contribute 
to the measured process. A precision lattice-QCD calculation can provide a cleaner gA directly 
from the SM and confront the experimental measurement, impacting many small-signal searches 
for SM-suppressed processes, as well as the astrophysical field for improving Solar-model physics. 

In early nucleon matrix element calculations, gA was considered the low-hanging fruit in lat- 
tice nucleon matrix-element calculations, the gold-plated equivalent of say f% or / = 2 %-% scatter- 
ing. The nucleon axial charge has been considered a benchmark quantity for lattice-QCD nucleon 
matrix elements. Due to limited resources, we were not able to provide sufficient statistics nor 
close-to-physical pion masses to better quantify our systematics; as a result, the overall uncertainty 
is quite large and overlaps the entire experimentally allowed range. However, many recent calcu- 
lations, involving multiple lattice spacings (i.e. continuum extrapolated) and improved statistics 
have been able to close the gap. Even as more computational resources become available, the best 
gA calculation is at 5% statistical error and a few standard deviations away from the experimental 
values. However, as seen on the right-hand side of Fig. 4, these data points are 10-15% lower 
than the experimental values, and those that seem to agree with experiment within 1-2 a tend to 
have large errorbars. A more reliable estimation of systematic uncertainties has not yet been made 
for the nucleon matrix elements. For example, estimates of the finite-volume effects might be too 
small. Before we jump to a conclusion concerning the existence of new physics or contemplate 
inherent problems with lattice QCD, we should reexamine all the systematics and push to control 
them at a reliable level of precision (say, percent-level or smaller). 

In the following, we consider the usual suspects for the systematics of a lattice measurement: 
renormalization, lattice discretization artifacts, excited-state contamination, chiral behavior and 
finite-volume corrections. 

4.1 Renormalization and lattice-discretization artifacts 

It is important to note that for the clover or 0(a) -improved Wilson fermions, we also need to 
0(a)-improve the axial current A^ = Za(cia +bAam q )(Ay L +acAd ll P) in the gA calculation, where 
the bA term vanishes in the chiral limit. Otherwise, the systematics due to lattice discretization arti- 
facts would come in at 0(Aqcd«) ~ 30% (for a = 0. 12 fm and Aqcd ~ 0.5 GeV). Mainz group has 
good control over the current and renormalization through Schrodinger-functional scheme, while 
LHPC does not. PNDME plans to make three lattice-spacing calculations to study and remove dis- 
cretization effects. The current-improvement 0(a) term, although suppressed for chiral fermions 
at finite a, should still be carefully considered as a potential residual chiral symmetry breaking; 
for example, in the case of domain-wall fermions, non-small m res can give rise to significant 0(a) 
artifacts. Note that when one composes ratios such as gA/ fn, it has been found that the ratio has 
better agreement with experiment at given lattice spacing. It is not only because the Za is canceled, 
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Figure 5: The bare nucleon charges as functions of source-sink separation (At € [0.96, 1 .44] fm) by PNDME 
(a « 0.12 fm, nijz ~ 310 MeV) with (9(4000) measurements. The black (blue) points are the charge values 
obtained from one-state (two-state) fits at each At, while the band indicates a simultaneous fit to all Af. 
Although the one- and two-state analyses give charges consistent within the statistical error, the two-state 
have smaller fluctuations about the central values and are more consistent with the simultaneous-fit analysis. 



but also the 0(a) systematics are reduced if there are remaining significant 0(a) artifacts in the 
currents. 

Among different lattice fermion actions, there are no observed significant artifacts on gA with 
lattice spacing below 0.12 fm. (See the right-hand side of Fig. 4 for lattice-spacing dependence. 
Data points are marked by different types of dashing for lattice spacings in the regions of 0.12, 0.09, 
0.06 fm.) It is safe to say that this is not the problem that causes the discrepancy from experiment 
(crisis) that we have observed, but it may become important and significant in the future when 
we try to become competitive with experiments to provide high-precision neutron-decay constants. 
More resources should be devoted to increasing the statistics, at least to the 0(1%) level, and we 
should not forget to revisit this systematic when higher-statistics measurements become available. 

4.2 Chiral extrapolation 

There is always a debate as to how reliable chiral effective theory (ChPT) is and how well 
it describes the data; there are diverse observations and nearly as diverse conclusions. Statistics 
matters when it comes to whether a lattice-QCD has successfully provide a reliable gA or not. In 
many of the chiral extrapolation shown on the left-hand side of Fig. 4, using either heavy-baryon 
chiral perturbation theory (HBChPT) or small-scale expansion (SSE) formulation, one needs to set 
a number of low-energy constants; those constants could be biased tuned to match the experimental 
values. We may need to consider using blinded techniques for extrapolation to avoid human bias 
in addition to making precise measurements and quoting systematics. Luckily, calculations on 
ensembles at the physical pion mass are an attainable near-term goal; however, they will require 
large-enough volumes for a gA calculation and long trajectories. 

4.3 Excited-state contamination 

Excited-state contamination poses a problem due to the nearby Roper resonance on top of the 
nucleon's poor signal-to-noise problem at large Euclidean time, which scales like g-Wv- 3 " 1 */ 2 )'. 
In terms of three-point correlators, from which we extract gA, to be confident that the signal is free 
of excited-state contamination, we need large source-sink separation; this large t also worsens the 



11 



Lattice Hadron Structure: Applications within and beyond QCD 



Huey-Wen Lin 



noise. The usual recipe is to check the source-sink dependence at low statistics and identify one 
value for full production; however, it is difficult to be completely sure that the chosen value of the 
separation would be clean of excited-state contamination at the full statistical precision. 

Strategies for dealing with the problem can be categorized in 3 main directions: First, even if 
we limit ourselves to single-state analysis, we can simply increase the source-sink distance (say to 
1.4 or 1.5 fm), where we can extract our favorite plateau with most of the kinematic and overlap 
factors canceling out, as done by RBC/UKQCD. Alternatively, we can include the excited-state 
degrees of freedom and extract the ground-state nucleon by fitting. This was first proposed and 
demonstrated in Ref. [31] in 2008. To facilitate accurate extraction of the excited state to disentan- 
gle it from the ground state, we can either add operators that overlap differently with the ground 
and excited states, or have high statistics to stabilize the analysis of multiple degrees of freedom. 
For example, see work done in Refs. [31, 11]. More groups have started to adopt similar ideas into 
their analyses, including LHPC, CLS/Mainz and ETMC. 

A recently proposed method called the "summation" method was introduced to take advantage 
of small-separation three -point correlators, extending the source-sink separation to where cleaner 
signal can be easily obtained given the same statistics. CLS/Mainz claim that the summation anal- 
ysis resolves the discrepancy in smaller-plateau single-state analysis; however, although Mainz's 
studies observe a noticeable upward shift in the central values, their final numbers are consistent 
with a single-plateau analysis at 1.1 fm or so. LHPC makes a similar observation; their final "sum- 
mation" analysis is consistent with the single-state analysis with largest separation. The method 
does not include both the transition matrix elements; when one works with nonzero momentum 
transfer, depending on the chosen kinematics, one should adopt a formulation that can pick up 
the dominant transition matrix and time factors so that excited-state transition form factors can be 
removed from the ground-state signal. 

It is worth mentioning that ETMC uses high statistics with multiple source-sink separations 
in the best study of excited-state contamination for nucleon matrix elements [32]. They increased 
the statistics at the large source-sink separation by an order of magnitude so that all the matrix el- 
ements have comparable statistical errors. ETMC also compared different analysis approaches, in- 
cluding summation (seven source-sink separations) and generalized eigenvector problem (GEVP), 
and found g& to have mild source-sink separation dependence as long as the nucleon operator is 
highly optimized to overlap with the ground-state; however, other matrix elements, such as the 
leading-moment of the quark momentum fraction, may not. Different quantities must be studied 
on a case-by-case basis. CSSM adopted the variational method on matrix-element correlators for 
the meson case; future work with nucleon matrix elements is in progress. 

To resolve the problem of excited-state contamination once and for all, the solution (in the au- 
thor's opinion) would be to fully include the excited-state degrees of freedom (including 2 nucleon- 
Roper transition matrix elements and even the Roper-Roper one) and improve the statistics to sta- 
bilize the analysis. An example is shown in Fig. 5 with preliminary PNDME measurements at 
0(4k) statistics; the 2-state analyses remain consistent even at smaller separation. One might be 
able to push to smaller separation (but it may depend on the smearing and operator choices); at 
some point, the second-excited state and yet higher states must come in. One still needs to proceed 
with caution and not too aggressively reduce the separation. RBC/UKQCD are adopting all mode 
averaging (AMA), and #QCD and JLQCD are using all-to-all approaches. Of course, an all-to-all 
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propagator would be very useful for constructing an excellent approximation to the nucleon ground 
state, as in variational methods. 

4.4 Finite-volume corrections 

The systematics caused by the finite-volume effects have been known to decrease the value of 
gA by significant amount. Two commonly used formulations to determine finite-volume corrections 
are HBChPT [33] and its variation SSE [20]. The finite- volume corrections in these references are 
very similar under the replacement of / = V2F K , ca = V2gAN, g\ = — £aa> near-physical A and 
g G A in the chiral limit. (A different formulation regarding finite-volume corrections by CSSM can 
be found in Ref. [34]). The left-hand side of Fig. 6 shows two different sets of parameter choices 
on multiple selected volumes L G [3,5] fm as function of pion mass. These parameters are set in 
Ref. [20] using a variety of constraints to F n = 86.2 MeV, ca = 1.5, gi = 2.6, g A = 1.15, while 
Ref. [33] use SU(6) relations to derive gA = 1 + (2/3) cos 2 y, gA/v = — 2cosi/f, #aa = —3. The 
finite-volume corrections at physical pion mass are 0-0.004 and 0.05-0.004 for volume changing 
from 3 to 5 fm (corresponding to m n L = 2.1-3.6), which would be well below the statistical error at 
the current level of lattice-calculation precision. Note that Ref. [20] has a larger volume correction 
mainly due to larger coupling in ca terms. One can also replace the formula with lattice m n /f„, but 
this only results in a small difference. Groups that estimate their finite-volume corrections using 
the Refs. [33, 20] have reported small systematics to less than 1% for m K L pa 3. 

However, an earlier quenched RBCK study [35] saw more significant effects, as shown in 
Refs. [33, 20]. Recent RBC/UKQCD Nf = 2 + 1 studies also see a significant central-value shift 
of about 0.05 and 0.12 with pion masses 670 and 420 MeV when the volume changes from 2.74 
to 1.82 fm (which corresponds to m n L going from 9.3 to 6.2 and 5.3 to 3.9, respectively). QCDSF 
has multiple lattice spacings {a £ [0.06,0.078] fm) with various pion masses and volumes, and also 
sees a bigger volume correction than suggested in Refs. [33, 20]. 

To compare the finite-volume corrections from ChPT and lattice data, we take the difference 
between two lattice gA values at fixed pion mass and lattice spacing but with differing volumes, and 
compare with the change predicted according to Ref. [20]; the right-hand side of Fig. 6 selects a 
subset of lattice calculation points from RBC/UKQCD and QCDSF. Note that although the lattice 
volume-shifts seems to be at the edges of their errorbars of the Ref. [20] line, they are systematically 
larger in magnitude with central values greater by a factor of 2 to 10. 

Finite- volume collections need to be done at smaller pion masses; many groups estimate the 
finite-volume correction at larger pion mass with the same m K L, but we have seen that yet larger 
m n L are needed for smaller pion masses. This effect could cause underestimation of the finite- 
volume corrections in gA- The volume changes can cause shifts of up to 0.3 in the lattice gA as 
suggested from selected lattice data in Fig. 6; this alone could explain the lattice and experiment 
discrepancy. We need to study these effects with greater statistics to rule out or better understand 
the form of the finite-volume effects. 

5. Summary and Outlook 

There have recently been many new developments, calculations of new quantities and new 
ideas for exploring and extending the reach of lattice QCD for hadronic structure. Many of the 
calculations presented here may not be final numbers from lattice QCD that can be quoted by 
experiments or model builders, and they may not be completely free from uncontrolled systematics. 
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Figure 6: (left) The finite-volume corrections for gA as a function of n-v\ with the parameter choices set 
in Refs. [33, 20] at volume set to 3, 4, 5 fm. Note that the ~\y has changed to n/3 to better match with 
experimental coupling values. The formulation suggests finite volume correction at physical pion mass for 
4 fm box is less than 1%. (right) The changes in gA as volumes changes at a fixed lattice spacing and 
various pion masses. The top panel selects examples of the QCDSF 2-flavor results, while the bottom one 
shows sample data from RBC/UKQCD. The lines indicate gA shifts corresponding to parameter choices 
in Ref. [20], which indicate larger finite-volume corrections. Note that although the lattice values agree 
within 1-2 sigma with Ref. [20], the central values are greater by up to an order of magnitude at the lightest 
QCDSF point or a factor of 2 for RBC/UKQCD. The volume corrections can cause up to 0.3 shift in these 
cases, which could account for the lattice-experimental discrepancy; we need to use high-statistics data to 
better understand this systematic. 

Even so, such preliminary explorations represent exciting progress, and we need people to make 
aggressive calculations that push the frontier. Hadronic structure is essential to understanding 
how the quarks and gluons make up hadrons. We push their precision so that we can make tests 
of the Standard Model when compared with reliable experimental results or evaluate the validity 
of nonperturbative QCD models. Hadronic low-energy couplings can be contributed to effective 
theories to extend the reach of theory to processes that are too complicated to directly calculate on 
the lattice, and they can be used in combination with precision low-energy experiments to probe 
new interactions induced by potential few-TeV BSM particles. We have moved from a period where 
we barely understood potential systematics to an era of examining them with greater statistics. By 
the application of increasing computational resources, we will push to larger volumes (M K L > 6) 
at physical pion mass and longer ensemble trajectories so that we will be able to compete with 
experiment in quantities such as g A . This is an achievable goal if we combine our efforts. 
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